First-order phase transitions in superconducting films: A Euclidean model 
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I. INTRODUCTION 

In the last few decades, a large amount of work has 
been done on the Ginzburg-Landau phenomenological 
approach to critical phenomena. An account on the state 
of the subject and related topics can be found, for in- 
stance, in Refs. USSSHllBEHIlS- Questions 
concerning the existence of phase transitions may also be 
raised if one considers the behavior of field theories as a 
function of spatial boundaries. The existence of phase 
transitions would be in this case associated to some spa- 
tial parameters describing the breaking of translational 
invariance, for instance, the distance L between planes 
bounding the system. Analyses of this type have been 
recently performed 0, 0, ^1 . In particular, if one con- 
siders the Ginzburg-Landau model confined between two 
parallel planes, which is assumed to describe a film of 
some material, the question of how the critical tempera- 
ture depends on the film thickness L can be raised. 

Studies on field theory applied to bounded systems 
have been done in the literature for a long time. In par- 
ticular, an analysis of the renormalization group in finite- 
size geometries can be found in 0,^3- These have been 
performed to take into account boundary effects on scal- 
ing laws. In another related topic of investigation, there 
are systems that present domain walls as defects, created 
for instance in the process of crystal growth by some pre- 
pared circumstances. At the level of effective field theo- 
ries, in many cases, this can be modeled by considering 
a Dirac fermionic field whose mass changes sign as it 
crosses the defect, meaning that the domain wall plays 
the role of a critical boundary separating two different 
states of the system 0, . Under the assumption that 
information about general features of the behavior of sys- 
tems undergoing phase transitions in absence of external 
influences (like magnetic fields) can be obtained in the 
approximation which neglects gauge field contributions 



in the Ginzburg-Landau model, investigations have been 
done with an approach different from the renormalization 
group analysis. The system confined between two paral- 
lel planes has been considered and using the formalism 
developed in Refs. [HII1II1, the way in which the criti- 
cal temperature is affected by the presence of boundaries 
has been investigated. In particular, a study has been 
done on how the critical temperature T c of a supercon- 
ducting film depends on its thickness L [TH . ITsl 191. In 
the present paper we perform a further step, by consid- 
ering in the same context an extended model, which be- 
sides the quartic field self-interaction, a sextic one is also 
present. It is well known that those interactions, taken 
together, lead to a renormalizable quantum field theory 
in three dimensions and which is supposed to describe 
first-order phase transitions. 

We consider, as in previous publications, that the sys- 
tem is a slab of a material of thickness L, the behavior 
of which in the critical region is to be derived from a 
quantum field theory calculation of the dependence of the 
renormalized mass parameter on L. We start from the 
effective potential, which is related to the renormalized 
mass through a renormalization condition. This condi- 
tion, however, reduces considerably the number of rele- 
vant Feynman diagrams contributing to the mass renor- 
malization, if one wishes to be restricted to first-order 
terms in both coupling constants. In fact, just two di- 
agrams need to be considered in this approximation: 
a tadpole graph with the 4> 4 coupling (1 loop) and a 
"shoestring" graph with the 4> 6 coupling (2 loops) (see 
Fig.l). No diagram with both couplings occur. The in- 
dependence appears from the treatment of the loop in- 
tegrals, as the material is confined between two planes 
a distance L apart from one another. We therefore take 
the space dimension orthogonal to the planes as finite, 
the other two being otherwise infinite. This dimension of 
finite extent is treated in the momentum space using the 
formalism of Ref . ^| • 
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The paper is organized as follows. In Section II, we 
present the model and the description of a bounded sys- 
tem through an adaptation of the Matsubara formalism 
suited for our purposes. The contributions from the 
two relevant Feynman diagrams to the effective poten- 
tial are established, as well as an expression showing the 
L-dependence of the critical temperature. In Section III, 
as we wish to compare our theoretical result with some 
experimental data, we need first to make a phenomeno- 
logical evaluation of the <f> 6 coupling constant, based on 
the analogous derivation made by Gorkov for the 4> A con- 
stant. The comparison with measurements is discussed 
in Section IV. Finally, in Section V we present our con- 
clusions. 



II. THE EFFECTIVE POTENTIAL IN THE cj> b 
GINZBURG— LANDAU MODEL 

We start by stating the Ginzburg-Landau Hamiltonian 
density in a Euclidean D-dimensional space, now includ- 
ing both tfi 4 and </> 6 interactions, in the absence of external 
fields, given by (in natural units, fi = c = ks = 1), 

H=\\V V \ 2 + \mlW\ 2 - + (1) 

where we are taking the approximation in which A > 
and r\ > are the renormalized quartic and sextic self- 
coupling constants. Near criticality, the bare mass is 
given by ttiq = a(T/To — 1), with a > and To being 
a parameter with the dimension of temperature. Recall 
that the critical temperature for a first-order transition 
described by the hamiltonian above is higher than To 
pof . This will be explicitly stated in Eq. J2U below. 
We consider the system confined between two parallel 
planes, normal to the distance L apart from 

one another and use Cartesian coordinates r = (x, z), 
where z is a (D — l)-dimensional vector, with correspond- 
ing momenta k = {k Xl q), q being a (D — l)-dimensional 
vector in momenta space. The generating functional of 
Schwinger functions is written in the form 

Z = J T><p*T>(pexp (- dx J d D - 1 zn(\ip\,\Vip\)^j , 

with the field ip(x,z) satisfying the condition of confine- 
ment along the x-axis, <p(x < 0, z) = <p(x > L, z) =const. 
Then the field should have a mixed series-integral Fourier 
representation of the form 

oo „ 

<p(x,m)= Cr. / «f Z, " 1 96(q)e- <w " x - <q "0(wn ) q), 

n— — oo 

(3) 

where oj n = 2?:n/L and the coefficients c n and 6(q) cor- 
respond respectively to the Fourier series representation 
over x and to the Fourier integral representation over 
the (D — l)-dimensional z-space. The above conditions 



of confinement of the x-dependence of the field to a seg- 
ment of length L allow us to proceed, with respect to 
the x-coordinate, in a manner analogous as is done in 
the imaginary-time Matsubara formalism in field theory 
and, accordingly, the Feynman rules should be modified 
following the prescription 

f dk x 1 ^ 2nir 

n— — oo 

We emphasize, however, that we are considering an Eu- 
clidean field theory in D purely spatial dimensions, so we 
are not working in the framework of finite-temperature 
field theory. Here, the temperature is introduced in the 
mass term of the Hamiltonian by means of the usual 
Ginzburg-Landau prescription. 

To continue, we use some one-loop results described in 
[TH IL?! l2l| , adapted to our present situation. These re- 
sults have been obtained by the concurrent use of dimen- 
sional and zeta-function analytic regularizations, to eval- 
uate formally the integral over the continuous momenta 
and the summation over the frequencies u) n . We get sums 
of polar (L-independent) terms plus L-dependent ana- 
lytic corrections. Renormalized quantities are obtained 
by subtraction of the divergent (polar) terms appearing 
in the quantities obtained by application of the modified 
Feynman rules and dimensional regularization formulas. 
These polar terms are proportional to T-functions hav- 
ing the dimension D in the argument and correspond to 
the introduction of counterterms in the original Hamil- 
tonian density. In order to have a coherent procedure in 
any dimension, those subtractions should be performed 
even for those values of the dimension D for which no 
poles are present. In these cases a finite renormalization 
is performed. 

In principle, the effective potential for systems with 
spontaneous symmetry breakin g is obtained, following 
the Coleman- Weinberg analysis 22] , as an expansion in 
the number of loops in Feynman diagrams. Accordingly, 
to the free propagator and to the no-loop (tree) diagrams 
for both couplings, radiative corrections are added, with 
increasing number of loops. Thus, at the 1-loop approxi- 
mation, we get the infinite series of 1-loop diagrams with 
all numbers of insertions of the 4 vertex (two external 
legs in each vertex), plus the infinite series of 1-loop di- 
agrams with all numbers of insertions of the </> 6 vertex 
(four external legs in each vertex), plus the infinite series 
of 1-loop diagrams with all kinds of mixed numbers of 
insertions of <f> and <f> 6 vertices. Analogously, we should 
include all those types of insertions in diagrams with 2 
loops, etc. However, instead of undertaking this compu- 
tation, in our approximation we restrict ourselves to the 
lowest terms in the loop expansion. We recall that the 
gap equation we are seeking is given by the renormaliza- 
tion condition in which the renormalized squared mass is 
defined as the second derivative of the effective potential 
U(ifo) with respect to the classical field (po, taken at zero 
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FIG. 1: Diagrams contributing to the renormalized mass at 
lowest order in the coupling constants. 



field, 



d 2 E%o) 



dtpo 



(5) 



<Po=0 



Within our approximation, we do not need to take into 
account the renormalization conditions for the interac- 
tion coupling constants, i.e., they may be considered as 
already renormalized when they are written in the Hamil- 
tonian. At the 1-loop approximation, the contribution 
of loops with only 4 vertices to the effective potential 
is obtained directly from as an adaptation of the 
Coleman-Weinberg expression after compactification in 
one dimension, 



D 



oo 

E 



(-l) s 



Si 2s 



2.s 



-9lPo 



iD-li 



(k 2 + an 2 + c 2 ) 



(6) 



In the above formula, in order to deal with dimension- 
less quantities in the regularization procedure, we have 
introduced parameters c 2 = m 2 / Air 2 /j, 2 , a — (L[i)~ 2 , 
<7i = — \/16n 2 u 4 ~ D and 4>o — Vol ^ D ^ 2 ^ 1 , where (po is 
the normalized vacuum expectation value of the field (the 
classical field) and fi is a mass scale. The parameter s 
counts the number of vertices on the loop. 

It is easily seen that only the s = 1 term contributes to 
the renormalization condition JSJ). It corresponds to the 
tadpole diagram. It is then also clear that all 6 -vertex 
and mixed </> - and <^> 6 -vertex insertions on the 1-loop di- 
agrams do not contribute when one computes the second 
derivative of similar expressions with respect to the field 
at zero field: only diagrams with two external legs should 
survive. This is impossible for a </> 6 -vertex insertion at 
the 1-loop approximation, therefore the first contribu- 
tion from the (f> 6 coupling must come from a higher-order 
term in the loop expansion. Two-loop diagrams with two 
external legs and only (j> A vertices are of second order in 
its coupling constant, and we neglect them, as well as all 
possible diagrams with vertices of mixed type. However, 
the 2-loop shoestring diagram, with only one </> 6 vertex 
and two external legs is a first-order (in 77) contribution to 
the effective potential, according to our renormalization 
criterion. 

Therefore the renormalized mass is defined at first or- 
der in both coupling constants, by the contributions of 



radiative corrections from only two diagrams: the tad- 
pole and the shoestring diagrams. The tadpole contribu- 
tion reads (putting s = 1 in Eq. |JBJ|), 



1 00 r 

U x {^L)=ii D ^- 9l 4>t E J 



iD-li 



(7) 



The integral on the D — 1 non-compactified momentum 
variables is performed using the dimensional regulariza- 
tion formula 



d d k 
k 2 + M ' 

for d = D — 1, we obtain 



r(i-|) 7t d / 2 _ 

M 1_d / 2 ' 



(8) 



[i D \/a- 



r (C-l)/2 



D 



E 



(an 2 +c 2)(3-z»)/2' 



(9) 



The sum in the above expression may be recognized as 
one of the Epstein-Hurwitz zeta-functions, Zf (^^;a), 



2j^ + l 4^-1 
2 2 7T 2 



which may be analytically continued to |2 

2^-3/2 



Z{ (u;a) 



OO 

+^E 



/ \ !- 2,y 








H) 


7 r 





' \ fjL 2 Ln 



1/2-1/ 



K u _ 1/2 (mnL) 



(10) 



where the K v are Bessel functions of the third kind. The 
tadpole part of the effective potential is then 



Ui(<p ,L) 



{2n) D/2 - 2 

OO 

+E 



D-2 / „ 



n=l 



m 



a 2 Ln 



D/2-1 



K D / 2 -i(mnL) 



(11) 



We now turn to the 2-loop shoestring diagram contri- 
bution to the effective potential, using again the Feyn- 
man rule prescription for the compactified dimension. It 
reads 

U 2 {^L) = (i 2D - 2 a^g 2 <P 2 

'(12) 

where g 2 = rj/l28ir 4 u 6 ~ 2D . After subtraction of the po- 
lar term coming from the first term of Eq. I|10|) we get 



f/r n) (^,i) 



4(2ir) D / 2 ^ \u 



E 



2 nL 



D/2-1 



K 



D/2-i(mnL) 
(13) 
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and 



rr(Ren), p _ Wo 

U 2 i<Po,L) - 



n=l 



\nLJ 



K D/ 2-x{mnL) 



occurs when all the three minima of the potential are 
simultaneously on the line U(ipo) = 0. This gives the 
condition 



(14) 



Thus the full rcnormalized effective potential is given by 

1 o t A a Tj k „mfin"l Tr (Ren) 



m {L) = w 



(18) 



t, t\ * 2 2 n 4 , V 6 i rr( Rcn ) 



(15) 



The renormalized mass with both contributions then 
satisfies an L-dependent generalized Dyson-Schwingcr 
equation, 



m 2 (L) = itiq — 



A 



2(2* 



E 



D/2-1 



Notice that the value m = is excluded in the above 
condition, for it corresponds to a second-order transition. 
For D = 3, which is the physically interesting situation 
of the system confined between two parallel planes em- 
bedded in a 3-dimensional Euclidean space, the Bessel 
functions entering in the above equations have an ex- 
plicit form, K 1 / 2 {z) — \f*e~ z /V2z, which replaced in 



l> i L^i \ n L) K D / 2 -ii mn L) Ea. (|16|) . performing the resulting sum, and reminding 



+ 



n 



4(2tt) d 



E 



KnLJ 



K 



D/2- 



_i(mnL) 



(16) 



Thus, the effective potential l|15f) is rewritten in the form 



that ml = a(T/T 



1), gives 



nr(L) = n[L-! 



AI ln(1 

8n L y 



+ 



TJTT 



8(2tt) 3 L 2 



ln(l 



,-m(L)L\ 



(19) 



E%o) 
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<Po, 



(17) 



where it is assumed that A, 77 > 0, a necessary condition 
for the existence of a first-order phase transition associ- 
ated to the potential (|17f) . Then, a first-order transition 



In Eq. (|18f) m(L) may have any strictly positive value and 
this condition ensures that we are on a point on the criti- 
cal line for a first-order phase transition. Then introduc- 
ing the value of the mass, Eq. l(TB|) , in Eq. (|T§|) , we obtain 
the critical temperature 



T C (L) = T c il - 1 



3A 2 V 
1677a / 



A 



8naL 



ln(l 



r 

16V L ) + 



>1 



64tt 2 <xL 2 



ln(l 



1 3\ 2 T 



(20) 



where 



T c = T 1 



3A 



, (21) 
1677a / 

is the bulk (L — > 00) critical temperature for the first- 
order phase transition. 



III. PHENOMENOLOGICAL EVALUATION OF 
THE CONSTANT 77 

Our aim in this section is to generalize Gorkov's 
microscopic derivation done for the \ip 4 model 
in order to include the additional interaction term r)(p 6 
in the free energy. Here, our interest is to determine the 
phenomenological constant 77 as a function of the micro- 
scopic parameters of the material, in an analogous way as 



it has been done for the constant A in the \<p 4 model. Us- 
ing Gorkov's equations combined with the self-consistent 
gap condition |25| the free energy density may be written 
in terms of the gap energy A(x) as 



/(A) 



N(0) 



£o|VA| 



T 



- 1 IA 



1674C(5K 4 
147 /7 4 4 C 2 (3) 



IAI 



(22) 

where N(0) is the density of states at the Fermi surface, 
£0 is the coherence length, vf the Fermi velocity and £(x) 
is the Riemann zeta-function. N(0) and £0 & r e given by 



2V(0) 



1 



47r 2 fcsTF V h J 



(PFY 



£0 « 0.13 



Tivf 



(23) 
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where Tp is the Fermi temperature and kg is Boltz- 
mann's constant. T is the temperature parameter intro- 
duced in Eq. (JJJ that can be obtained from the first-order 
bulk critical temperature by means of Eq. (|21[1 . Introduc- 



In order to be able to compare our results with some 
experimental observation, we should restore SI units (re- 
member that so far we have used natural units, c — h — 
kg = 1) |2^. In SI units, the exponent in the partition 
function © has a factor 1/fcflTb. Then, we must divide 
by kgT the free energy density in Eq. l|24[l . Moreover, we 
rescale the fields and coordinates by ip ncw — W^o/kFsTQip 
and x ncw = x/£oi which gives the dimensionless energy 
density and, comparing with Eq. (JTJ, we can identify the 
phenomenological dimensionless constants A, 77, and mo, 
with a = 1 [26J, 

A « 111.08 (^j , r) « 0.04257 A 2 , m 2 = £ - 1. 

° (25) 

By replacing the above constants in Eq. J5DJ), we get 
the critical temperature as a function of the film thick- 
ness and in terms of microscopic tabulated parameters 
for specific materials. 

IV. COMPARISON WITH EXPERIMENTAL 
DATA 

We remark that Gorkov's original derivation of the 
phenomenological constants is valid only for perfect crys- 



« t '= T -{ 1 -( 1 + S) 



We consider that other effects, such that of the sub- 
strate over which the superconductor film is deposited, 
should be taken into account. In the context of our 
model, however, we are not able to describe such effects 
at a microscopic level. We therefore assume that they 



ing the order parameter tp = y / 2~/V(0)^oA in Eq. 122(1 we 
obtain 



(24) 



I 

tals, where the electron mean free path I is infinite. 
However, we know that in many superconductors the 
attractive interaction between electrons (necessary for 
pairing) is brought about indirectly by the interaction 
between the electrons and the vibrating crystal lattice 
(the phonons). Considering that this interaction will be 
greater if we have impurities within the crystal lattice, 
consequently the electron mean free path is actually fi- 
nite. The Ginzburg-Landau phenomenological constants 
A and r\ and the coherence length are somehow related to 
the interaction of the electron pairs with the crystal lat- 
tice and the impurities. A way of taking these facts into 
account preserving the form of the Ginzburg-Landau free 
energy is to modify the intrinsic coherence length and 
the coupling constants. Accordingly |2q|. £0 - * J " 1 '' 2 £o> 
A — > 2r~ 3 / 2 A and rj — > 4r~ 3 ?7, where r ps 0.18i? _1 , with 
R = to/l. Then, Eq. (HU becomes 



will be translated in changes on the values of the cou- 
pling constants A and 77. So, we propose as an Ansatz 
the rescaling of the constants in the form A — ► aX and 
rj — > a 2 rj. We may still combine both parameters R and 
a as r = aR. Eq. l]2(:iH is then written as 
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4n'v F e N(o)\ 



1674 C(5) 
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T C (L) = T c {l - 1 
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(27) 
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FIG. 2: Critical temperature T C (K) as function of thickness 
L (A), from Eg. l27H and data from Ref. [27J for a supercon- 
ducting film made from aluminum. 



In Fig. 2 we plot Eq. (|27[l to show the behavior of 
the transition temperature as a function of the thickness 
for a film made from aluminum. The values for Al of the 
Fermi temperature and the bulk critical temperature are 
T F = 13.53 x 10 4 K and T c = 1.2 K, respectively. We 
see from the figure that the critical temperature grows 
from zero at a nonnull minimal allowed film thickness 
above the bulk transition temperature T c as the thickness 
is enlarged, reaching a maximum and afterwards start- 
ing to decrease, going asymptotically to T c as L — > oo. 
We also plot for comparison some experimental data ob- 
tained from Ref. j2?|]. We see that our theoretical curve 
is in qualitatively good agreement with measurements, 
especially for thin films. 

The experimental evidence showing that in some su- 
perconducting films the transition temperature is well 
above the bulk one has been reported in the literature 
since the 1950s and 60s [H El H EH- On the the- 
oretical side, a formula for the transition temperature 
was written within BCS theory in terms of the electron- 
phonon coupling constant, the Debye temperature and 
the Coulomb coupling constant 32] . This formula was 
used to explain observed increases in the critical tem- 
perature of thin composite films consisting of alternating 
layers of dissimilar metals |2Jj • In Ref. [33] a molecular- 
dynamic technique was applied to obtain the phonon fre- 
quency spectrum which led to the same results. Mech- 
anisms accounting for the sharp drop in T C (L) for very 



thin films were also discussed in Ref. |27|. The authors 
conclude that the most important influence on T C {L) was 
the interaction of the film with the substrate, described 
by a model in Ref. |34| . 

It is interesting that in recent reports on copper oxide 
high-transition-temperature superconductors, the criti- 
cal temperature depends on the number of layers of Cu02 
in a similar way as above: first it rises with the number 
of layers and, after reaching a maximum value, then de- 
clines. See [35| and references therein. 

This behavior may be contrasted with the one shown 
by the critical temperature for a second-order transition. 
In this case, the critical temperature increases monoton- 
ically from zero, again corresponding to a finite minimal 
film thickness, going asymptotically to the bulk transi- 
tion temperature as L — > oo. This is illustrated in Fig. 
3, adapted from Ref. |36|. with experimental data from 
|37|. (Such behavior has also been experimentally found 
by some other groups for a variety of transition-metal ma- 
terials, see Refs. E3,E3,E3-) ^ nce * n ^ e P reseil t wor k a 
first-order transition is explicitly assumed, it is tempting 
to infer that the transition described in the experiments 
of Ref. (23 is first order. In other words, one could say 
that an experimentally observed behavior of the critical 
temperature as a function of the film thickness may serve 
as a possible criterion to decide about the order of the 
superconductivity transition: a monotonically increasing 
critical temperature as L grows would indicate that the 
system undergoes a second-order transition, whereas if 
the critical temperature presents a maximum for a value 
of L larger than the minimal allowed one, this would be 
signalling the occurrence of a first-order transition. 



V. CONCLUSIONS 

As seen in previous works, a superconducting system 
confined in some region of space may lose its character- 
istics if the dimensions of this region become sufficiently 
small. This is due to the fact that the critical tempera- 
ture depends on these dimensions in such a way that it 
vanishes below some finite minimal size. This has been 
verified in a field-theoretical framework for a Ginzburg- 
Landau model describing a second-order phase transi- 
tion. In the present paper, we have studied the criti- 
cal temperature behavior of a sample of superconducting 
material in the form of a film, but we have included in 
the model a (f> 6 self-interaction term, thus implying that 
we are now dealing with a first-order transition. In the 
case we have treated, a sharply contrasting behavior of 
the critical temperature, as a function of the film thick- 
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Tc vs L for superconducting Nb film 

9-5 -i Ref. [26] 

• Data from Ref. [27] 
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FIG. 3: Critical temperature T c (K) as a function of the thick- 
ness L (A) for a second-order transition, as theoretically pre- 
dicted in Ref. [3(J. Dots are experimental data taken from 
Ref. [37| for a superconducting film made from niobium. 



ness, was obtained with respect to the corresponding one 
for a second-order transition. This possibly indicates a 
way of discerning the order of a superconducting transi- 
tion from experimental data, according to the profile of 
the curve T c vs L. 

Also importantly, for our derivation of the first- 
order transition critical temperature curve, we needed to 
phenomenologically evaluate the 6 coupling constant, 
which, as far as we know, is not present in the literature. 

Finally, we also remark that in D = 3, for second-order 
transitions, one considers m — and that leads to the 
need of a pole-subtraction procedure for the mass |10| . 
In our case such a procedure is not necessary, as a first- 
order transition must occur for a non-zero value of the 
mass. This fact, together with the closed formula for the 
Bessel function for D = 3, allows us to obtain the exact 
expression (|20H for the critical temperature. 



Acknowledgments 

This work has received partial financial support from 
CNPq and Pronex. 



L. Halperin, T. C. Lubensky, S.-K. Ma, Phys. Rev. Lett. [21 
32, 292 (1974). 

I. Affleck and E. Brezin, Nucl. Phys. 257, 451 (1985). [22 
I.D. Lawrie, Phys. Rev. B 50, 9456 (1994). 

I.D. Lawrie, Phys. Rev. Lett. 79, 131 (1997). [23 

E. Brezin, D.R. Nelson and A. Thiaville, Phys. Rev. B 

31, 7124 (1985). [24 

L. Radzihovsky, Phys. Rev. Lett. 74, 4722 (1995). [25 

M.A. Moore, T.J. Newman, A.J. Bray and S.-K. Chin 

C. de Calan, A.P.C. Malbouisson and F.S. Nogueira, 

Phys. Rev. B 64, 212502 (2001). [26 

A.P.C. Malbouisson, F.S. Nogueira and N.F. Svaiter, 

Eur. Phys. Lett. 41, 547 (1998). 

L.M. Abreu, A.P.C. Malbouisson, I. Roditi, Physica A [27 
331, 99 (2004). 

A.P.C. Malbouisson and J.M.C. Malbouisson, J. Phys. [28 
A: Math. Gen. 35, 2263 (2002). 

A.P.C. Malbouisson, J.M.C. Malbouisson and A.E. San- 
tana, Nucl. Phys. B 631, 83 (2002). [29 
A.P.C. Malbouisson, J.M.C. Malbouisson, A.E. Santana 
and F.C. Khanna, Mod. Phys. Lett. A 20, 965 (2005). 
J. Zinn- Justin, Quantum Field Theory and Critical Phe- [30 
nomena (Clarendon Press, Oxford, 1996), chapter 36. 
J.L. Cardy (ed.), Finite-Size Scaling (North Holland, [31 
Amsterdam, 1988). 

CD. Fosco and A. Lopez, Nucl. Phys. B 538, 685 (1999). [32 
L. Da Rold, CD. Fosco and A.P.C. Malbouisson, Nucl. [33 
Phys. B 624, 485 (2002). 

L.M. Abreu, A.P.C. Malbouisson, J.M.C. Malbouisson, [34 
A.E. Santana, Phys. Rev. B 67, 212502 (2003). 
A.P.C. Malbouisson, Phys. Rev. B 66, 092502 (2002). [35 
M. Le Bellac, Quantum and Statistical Field Theory, Ox- 
ford University Press, Oxford, 1991. 



G. N.J. Ananos, A.P.C. Malbouisson and N.F. Svaiter, 
Nucl. Phys. B 547, 221 (1999). 

S. Coleman and E. Weinberg, Phys. Rev. D 7, 1888 
(1973). 

A. Elizalde and E. Romeo, J. Math. Phys. 30, 1133 
(1989). 

LP. Gorkov, Sov. Phys. JETP 9, 636 (1959). 

A. A. Abrikosov, L. P. Gorkov, and I. E. Dzyaloshinskii, 

Methods of Quantum Field Theory in Statistical Physics, 

1963. 

H. Kleinert, Gauge Fields in Condensed Matter, Vol. 1: 
Superflow and Vortex Lines, World Scientific, Singapore, 
1989. 

M. Strongin, R.S. Thompson, O.F. Kammerer and J.E. 
Crow, Phys. Rev. B 1, 1078 (1970). 

V. Buckel and R. Hilsch, Z. Phys. 132, 420 (1952); 138, 
109 (1954); N.V. Zavaritsky, Dokl. Akad. Nauk SSSR 86, 
501 (1952). 

M. Strongin, O.F. Kammerer, J.E. Crow, R.D. Parks, 
D.H. Douglass, Jr. and M. A. Jensen, Phys. Rev. Lett. 
21, 1320 (1968). 

M. Strongin and O.F. Kammerer, J. Appl. Phys. 39, 2509 
(1968). 

M. Strongin, O.F. Kammerer, D.H. Douglass, Jr. and 

M.H. Cohen, Phys. Rev. Lett. 19 121 (1967). 

W.L. McMillan, Phys. Rev. 167, 331 (1968). 

J.M. Dickey and A. Paskin, Phys. Rev. Lett. 21, 1441 

(1968). 

L.N. Cooper, Phys. Rev. Lett. 6, 689 (1961); W. Silvert 

and L.N. Cooper, Phys. Rev. 141, 336 (1966). 

M.V. Ramallo, Eur. Phys. Lett. 65, 249 (2004); S. 

Chakravarty, H.-Y. Kee and K. Volker, Nature 428, 53 

(2004). 



8 



[36] L.M. Abreu, A.P.C. Malbouisson and I. Roditi, Eur. 

Phys. J. B 37, 515 (2004). 
[37] J. Kodama, M. Itoh and H. Hirai, J. Appl. Phys., 54, 

4050 (1983). 

[38] H. Raffy, R.B. Laibowitz, P. Chaudhari and S. Maekawa, 
Phys. Rev. B 28, 6607 (1983). 



[39] M.S.M. Minhaj, S. Meepagala, J.T. Chen and L.E. 

Wenger, Phys. Rev. B 49, 15235 (1994). 
[40] A.V. Pogrebnyakov, J.M. Redwing, J.E. Jones, X.X. Xi, 

S.Y. Xu, Qi Li, V. Vaithyanathan and D.G. Schlom, 

Appl. Phys. Lett. 82 (2003) 4319. 



